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We discuss mechanical buckling instabilities of a rigid film under compression interacting repul- 
sively with a substrate through a thin fluid layer. The buckling occurs at a characteristic wavelength 
that increases as the l/4th power of the bending stiffness, like a gravitational instability studied 
previously by Milner et al. However, the potential can affect the characteristic buckling wavelength 
strongly, as predicted by Huang and Suo. If the potential changes sufficiently sharply with thick- 
ness, this instability is continuous, with an amplitude varying as the square root of overpressure. 
We discuss three forms of interaction important for the case of Langmuir monolayers transferred to 
a substrate: Casimir-van der Waals interaction, screened charged double-layer interaction and the 
Sharma potential. We verify these predictions numerically in the Van der Waals case. 

PACS numbers: 46.32,+x, 46.70, 64.70 



I. INTRODUCTION 

The advent of controlled molecular scale films and de- 
position methods has revealed a number of fine-scale 
wrinkling instabilities dii. At the same time, several 
new and general features in the buckling of macroscopic 
films have been identified [3,|5, 6]. Some of these are elab- 
orations of the simple Euler buckling of a compressed rod 
or sheet [3, IE] • Recently Cerda and Pocivavsek Q have 
considered rigid, compressed sheets on the surface of a 
liquid in the presence of gravity. This extends an earlier 
treatment of Milner, Joanny and Pincus 10] adapted to 
lipid monolayers at an air-water interface. Under these 
conditions the Euler buckling occurs not at zero wavevec- 
tor but at a finite wavevector determined by the bending 
modulus and liquid density. We call this mode of buck- 
ling gravity-bending buckling. 

Folding structure of laterally-compressed surfactant 
monolayers at the air-water interface is a well-known phe- 
nomenon [ill, E3] • The initial instability leading to these 
folds may be related to the gravity-bending buckling 
noted above. The observed folding length scale resem- 
bles the predicted wavelength of gravity-bending buck- 
ling . Analogous folding has recently been observed in 
solid nanocrystal monolayers EH- 

Additional topographic structure is observed when 
compressed lipid monolayers are transferred to a solid 
substrate via the inverted Langmuir-Schaeffer method 
[bH ]. These supported monolayers and bilayers are in- 
creasingly common in the study of biological membranes 
EE EE 53, EE El- The layer thus transferred is posi- 
tioned for easier study. Initially these transferred lay- 
ers are separated from the substrate by a cushion of 
the carrier liquid. Any topographic patterning of these 
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transferred layers can readily be observed [lg, [20( . Such 
patterning is likely affected by interaction with the sub- 
strate. Likewise, any buckling of a supported monolayer 
must be affected by the substrate. To study transferred 
monolayers under the high compressions where buckling 
is expected seems feasible, though to our knowledge no 
such studies have been performed. 

In this paper we investigate a class of wrinkling insta- 
bilities that are generalizations of gravity-bending buck- 
ling. These instabilities occur when a deformable surface 
such as a lipid monolayer lies above a solid substrate on 
a cushion of fluid. The interaction of the surface with 
the substrate can then play the role of gravity. This 
interaction alters the gravity-bending instability in sev- 
eral ways. It makes the unstable wavelengths depend on 
depth d. These wavelengths generally are much smaller 
than those predicted by the gravity-bending instability. 

In order to investigate the effect of the substrate in its 
simplest form, our treatment neglects several effects that 
may be important in practice [2 11 ] - In practice some ex- 
ternal forcing on the film is required in order to create the 
lateral pressure to induce buckling. For example regions 
outside the region under study may be bonded to the sub- 
strate. In practice the time dependence of the buckling 
may be important in determining its wavelength. This 
is especially true in cases where the equilibrium buck- 
ling transition is discontinuous. These effects have been 
extensively explored in the semiconductor film literature 
[22l I23I Hj] . Below we shall merely assume that a uni- 
form unixial stress is imposed on the region in question 
and will only investigate the initial buckling instability 
as influenced by the substrate interaction. Further, we 
shall suppose that the buckling film is inextensible. In 
the appendix we show that the inextensible approxima- 
tion is appropriate for the lipid monolayers such as those 
of Ref. 11 and 12. 

Moreover, the interaction alters the qualitative na- 
ture of the instability. The gravity-bending instability 
is a runaway or subcritical instability at constant surface 
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pressure. Here the amplitude of the wrinkles jumps from 
zero to a large value determined by other aspects of the 
system. However, substrate interactions can change this 
behavior, making the amplitude a continuous function of 
surface pressure. The criterion for a continuous transi- 
tion can be stated generally in terms of the second and 
fourth derivative of the interaction potential with sepa- 
ration. 

In the following sections we discuss the substrate- 
bending instability in terms of a general interaction po- 
tential 4>{d). We determine the unstable wavelength and 
the condition for a continuous transition. In the following 
section we consider three specific potentials commonly 
encountered in liquid films. The first is the Casimir-van 
der Waals interaction: <fr(d) — A/ (\2ird 2 ). The second 
is the screened charged double-layer interaction common 
in aqueous films with charged interfaces. The third is 
the Sharma potential used to describe molecularly thin 
water films. All of these potentials produce continuous 
wrinkling for sufficiently thick films. We conclude that 
this substrate-bending buckling should be readily observ- 
able. 
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(a). Flat state 



(b). Buckled state 



FIG. 1: One dimensional model used in the theory. The y- 
axis is pointing out of the figure, (a). Initial Flat state with 
no compression, (b). Buckled state with large enough com- 
pression. 



II. THEORY OF WRINKLING INSTABILITY 

A. Wavelength of Microscopic Wrinkles 

In the analysis below, we will consider a simplified 
model of a Langmuir monolayer. Suppose an insoluble 
surfactant layer is sitting at the interface between air 
and a liquid subphase. The elastic property of a layer, 
with finite thickness t, is characterized by the bending 
modulus B HI]: 



and relax the strain energy in a third dimension, similar 
to the Euler buckling of a rod. We recall that the crit- 
ical pressure n c is independent of the compressibility of 
the layer, although the corresponding strain depends on 
the compressibility 7]. In the following analysis a con- 
stant external pressure H ex is exerted on the Langmuir 
monolayer. A buckling transition is induced by displac- 
ing the boundary. In the wrinkled state, total area S of 
the surfactant layer is: 
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where E is the Young's modulus and v is the Poisson ra- 
tio in the continuum theory. A solid substrate is placed 
under the subphase, as shown in Fig[TJ Interaction be- 
tween the solid substrate and liquid subphase takes the 
form of a substrate potential </>(d) [3l|. For different in- 
teractions, the functional forms of (f>(d) are different [32| . 
For example, if the interaction is of pure van der Waals 
type, we have: 
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where A is the Hamaker constant [3l|. It can take pos- 
itive or negative values depending on properties of sub- 
strate and subphase. Possible retardation effects are not 
considered in this paper. In a Langmuir trough, one can 
compress the surfactant layer with external pressure H ex . 
Elastic strain energy is stored in the elastic layer upon 
compression. It is expected that if H ex exceeds some crit- 
ical value n c , the elastic layer will enter a buckled state 
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The above integral is taken over the projected area So 
on the horizontal x-y plane and ds is a surface element 
in So- The quantity £,(x,y) is the vertical displacement 
of interface from a flat state. The height profile of the 
interface is: d{x, y) = do + £,{x, y), where do is the height 
of a flat state with no surface deformation. We assume 
that no subphase fluid enters or leaves, so that the volume 
of subphase under the initial flat surface is fixed during 
deformation: 



f £(x,y)ds = 0. 

J So 



(4) 



Another parameter needed is the surface density of sur- 
factant: a = N/S, where N is the total number of sur- 
factant molecules. For constant external pressure H ex , 
the Gibbs free energy is written as (nj HH : 



G = 7o (S - S ) +F l +F b + F l + U ex S ■ 



(5) 
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where 70 is the surface tension of a free interface without 
any compression. The first term denotes the change in 
interfacial energy. F\ is the surfactant free energy. It is 
related to a via the relation: Fi = Sfi(a). The last term 
is an analogy of the pressure P times volume V term 
in the Gibbs free energy of a conventional gas. Fb and 
Fi are bending energy and substrate potential energy, 
respectively. For small surface deformation (|£/do| <C 1), 
the bending energy Fb is given by the Helfrich energy 

MM- 



In the above expression, £ denotes the derivative of £ 
with respect to the x-coordinate. In the one dimensional 
model, the mean curvature C is given by: 

c o (1 • i 2 ) " ' • (10) 

Using this expression the bending free energy(per unit 
length in y-direction) fb and its expansion in small de- 
formation approximation take the form: 
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(C-C ) 2 (1 + (V£) 2 ) 1/2 ds, (6) 



where C(x,y) is the mean curvature of the interface 
and Co is the spontaneous curvature of surfactant layer. 
Again we have expressed the integral in terms of the pro- 
jected area as in equation ([3]). Co has no effect on the 
following analysis and we shall neglect it henceforth [28| . 
The substrate potential energy is: 



Fi 
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where the initial flat state is chosen as reference state for 
the potential energy. Our treatment can be simplified 
by supposing that the film is inextensible, as justified in 
the appendix. In that case, the molecular density a is 
fixed and the energy F\ is a mere constant. Then for 
different values of external pressure Tl ex , the equilibrium 
configuration of the system is obtained by minimizing the 
Gibbs free energy with the inextensibility constraint of 
the surfactant layer: S — const. Introducing a Lagrange 
multiplier 9 to incorporate this constraint, the functional 
that we need to minimize is: 



G' = G-9 \S 
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In the rest of this paper, we assume that relaxation 
of strain energy only occurs in the x-direction, as shown 
in Figfl] In the y-direction, the system has translational 
invariance. As a result, we may minimize the Gibbs free 
energy per unit length in y-direction, which we denote as 



/ = g-eU- J l (i + m^dx^j 

= 7o (L - L ) + Lfi(a) + f b + fi + U ex L 



1- 1 a + dry^dx 

Lo 



(9) 



where L is the total length of surfactant layer and Lo 
is the projected length in x-direction, as shown in FigfT] 
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Similarly, the substrate potential energy /j and its ex- 
pansion take the form: 
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4>(d )L. (12) 



Here (f>'(do), 4>^ (do) . . . denote the derivatives of <f> with 
respect to z-coordinate and are evaluated in the flat state. 
Integration of the first order term in £ vanishes because 
of volume conservation, equation Furthermore, we 
will choose the origin of x-coordinate such that integra- 
tion ranges from —Lo/2 to Lq/2. Different coordinate 
systems differ only in small boundary terms, which are 
negligible if the system is large enough. In the following 
discussion it will be clear in what sense we mean by large 
enough. Minimizing g' with respect to Lo, L and surface 
undulation and keeping only the lowest order terms 
in £ , we get the following equilibrium equations of state 



B 



d 4 e 

dx 1 



Hex ~ 70 + <£(do) + 6 

70 - (f>{do) + v H 777TT — 

d(L) 

(0 + 0(do))£ + ^ (2) (doK 



= 
= 
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(13) 
(14) 

(15) 



The Lagrange multiplier 6 is related to the external pres- 
sure H ex '- 9 = 70 — Hex — 4>(do)- Moreover, the equilib- 
rium shape of the interface £(x) must satisfy the above 
differential equation. Using an Ansatz of sinusoidal de- 
formation: £(x) = hsbx(qx), where h is the amplitude 
and q is the wavenumber, we obtain: 
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Bq 4 + { lQ -Il ex )q 2 + <f>( 2 \d ) = 0, (16) 

where 9 = 70 — Tl ex — (f>(do) has been used (by equation 
(fl~3|) ). We have the following relation between external 
pressure H ex and wavenumber q: 



classical theory of phase transition [30( , a first-order (dis- 
continuous) transition will occur if the coefficient of the 
fourth order term of free energy expansion with respect 
to the order parameter (£(x) in our case) is negative. In 
this section, we will include terms up to the fourth or- 
der of £. The inextensibility constraint of the surfactant 
layer yields: 



70 = Bq 2 
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An equivalent equation was obtained by Huang et al [2!| • 
Minimizing the right hand side of equation (|17[) with re- 
spect to q, we get the smallest external pressure n c for 
buckling instability and critical wavenumber q c : 



j \ 1/2 
£ 2 ) dx 



const . 



(20) 



Under the inextensibility constraint (|20p , we can drop 
constant terms in expression ([9]) and minimize with re- 
spect to the following functional of Lq and £(x): 
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Ac 

n, 



D 



1/4 



,9i = (U ex - 7o)i + h + fi + <l>{do)L . 
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If the gravitational energy of a liquid subphase (in air) 
is considered, we may choose the flat state as reference 
and take (j>(d) = pg£ 2 /2. The threshold external pressure 
II C and critical wavenumber q c in this particular case are: 



n, 
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p 5 \i/4 
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The results in equation (fl9|) were obtained by Milner 
et al. in the paper [Icj . Gravitational energy is impor- 
tant in macroscopic scale with relatively thick liquid sub- 
phase. However, in microscopic scale (d <100nm), differ- 
ent types of substrate interaction between solid substrate 
and liquid subphase, e.g. van der Waals interaction, be- 
come dominant, while gravity is negligible. The buckling 
transition now leads to microscopic wrinkles. Our result 
in equation (|18|) is a generalization to this microscopic 
range. The effect of different types of interaction will be 
discussed in a later section. 



B. Second-order Buckling Transition 

In the previous section, we used the small deforma- 
tion approximation and expanded the Gibbs free energy 
to the lowest order in surface displacement £(x). A re- 
lation between external pressure Tl ex and wavenumber 
q of wrinkles is obtained in equation (|17[) . In order to 
study the undulation amplitude h and possible order of 
the buckling transition, higher order terms in the expan- 
sion should be included in the analysis. We are partic- 
ularly interested in finding out the existence conditions 
for a continuous, second-order transition. From Landau's 



The projected length Lq in x-direction and the surface 
deformation £(x) are not independent variables. They 
are related through the constraint ([2H)l . Assuming si- 
nusoidal deformation: £(x) = hsin(qx), it leads to an 
expression of Lq in terms of L, h and q. We will revisit 
the assumption of incompressibility in the Discussion sec- 
tion. Inserting this expression for Lq into equation (|2ip . 
we see that the functional g\ has the form of Landau 
free-energy expansion [3fj| . Firstly, we compute the ex- 
pression for Lq. Expanded to the fourth order in £(x), 
i.e. hsm(qx), the constraint (fJD]) is: 
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In the last step, we keep only extensive terms, propor- 
tional to the size of the system Lq, and neglect boundary 
terms. The boundary terms are of the order of a wave- 
length of wrinkles: 2n/q. The approximation made in 
equation (f2"2"|) is thus essentially that the wavelength of 
wrinkles is much smaller than the dimension of Langmuir 
system in the x-direction: \qL\ ^> 1. If this condition is 
satisfied, the approximation (|22[) is valid and we get the 
following expression for the projected length Lq: 
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where the slope amplitude h — hq is a dimensionless pa- 
rameter and \ h\ <C 1 in the small deformation approxima- 
tion. Similarly, expanding the functional g\ in equation 
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(121) to the fourth order in £ and neglecting boundary 
terms, we get: 



91 



h 2 , . 



<Xdo 



<^> (2) W 



64 



lOBg 2 - 30(d o ) 



(24) 



where = 70 — IT^ — 0(<io) is the Lagrange multiplier 
defined in the last section. g\ depends on the dimension- 
less variable h and wavenumber q. Using the expression 
of Lq in equation (|23|) : 
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The equilibrium configuration minimizes the value of 
g 2 . In the above expression g 2 — > — 00 as q — > 00. It seems 
that the minimum value of 32 doesn't exist. This paradox 
is solved by noticing that the value of wavenumber q can't 
vary arbitrarily. In the small deformation approximation, 
we have assumed that h = hq is a small quantity. As a 
result, q can't be arbitrarily large. Furthermore, q should 
be close to its critical value q c near transition: 
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where 8 is a small dimensionless parameter and 8 > if 
wavenumber q is a continuous function of compression. 
Using the expression (|26p for q and keeping the terms up 
to the second order in S, we obtain: 
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The minimum value of 32 exists for some positive x 
and 8, only if the coefficients A\ and A2 satisfy the in- 
equalities: 



A 1 < ; A 2 > . 



(31) 



Furthermore, if inequalities (|3ip are true, gi can achieve 
its minimum at h m i n and S m i n : 
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where D~ h and D$ are some positive coefficients. Obvi- 
ously, h m i n and S m i n are positive when n ea; > II C and 
approach zero, as the external pressure Ii ex approaches 
its threshold value Il c from above. In other words, the 
values of h m i n and S m i n can be made arbitrarily small 
by approaching transition point. As a result, the terms 
involving 8 in A\ and A2 are of higher order and can be 
neglected in comparison with the finite constant term. 
The inequalities (|3"Tj) are reduced to: 



(33) 
(34) 



n c -n ex < 

a - 166 - 70 - 70(do) > 0. 



The first inequality will be true if Tl ex > II C . Using the 
expression of 7, a and 6, we can rewrite the second in- 
equality as: 



Be 



(d )) /^ {2) (d )-2 (b^ (d )) 1/2 > 7(n c -n e:c ). 

Since the right hand side approaches zero from below as 
Ii ex decreases to II C , we get the condition: 
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where the coefficients A\ and A 2 are quadratic functions 
of 8: 
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where the two parameters 6 and a are defined as: 
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where we have used the expression (jT8|) for q c . In 
order to have a second-order buckling transition, the 
inequality ()36|) is the condition that must be satisfied 
by the substrate potential <j)(d). The above criterion 
for the order of buckling transition doesn't change 
after including the second harmonic term ah $m(2qx). 
Because a doesn't affect A\ and enters A 2 as a 2 h 2 , the 
second harmonic term affects g 2 only at order h 6 or 
higher. As a result, it does not change the criterion in 
inequality (|36p . In the case of gravitational potential 
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energy of liquid subphase: <f){d) — pg£ 2 /2, we have 
(j) {2) (d ) = pg > and <p {4) {d ) = 0. The relation §jc$ 
cannot be satisfied. Therefore, we will have a first-order 
buckling transition in this case [ic| |. However, the 
relation (j36|) may be satisfied for some types of substrate 
potentials. We will discuss its explicit forms in the next 
section and consider the possibility of a second-order 
buckling transition for different types of interaction. 
Although our conditions are necessary for a continuous 
transition, they are not completely sufficient. To show 
that no discontinuous transition occurs, we would have 
to show that no displacement £(2;) has a H m i n < n c . 
We considered only small and harmonic £(x). 



liquid subphase do, van der Waals potential energy den- 
sity decreases as Ajd\. Meanwhile, the local curvature 
of surfactant layer C is of the order of 1/do, so bending 
free energy density varies as B/dq, which has the same 
functional form as van der Waals interaction. As a result, 
we expect the criterion doesn't depend on the thickness 
of the liquid subphase do and is a relation between the 
Hamaker constant A and bending modulus B. As long 
as A and B satisfy condition (|39|) . microscopic wrinkles 
can be formed through a second-order buckling transi- 
tion. For a lipid monolayer B > lOfcT 36], which is 
4 x 10~ 20 J at 25 centigrade. The smallest A c compatible 
with equation ([39]) is: 



III. EXAMPLE POTENTIALS 

In this section, we will consider some examples of sub- 
strate potentials and find out the corresponding existence 
conditions for second-order buckling transition. 

A. Non-retarded Van der Waals Interaction 

Non-retarded van der Waals interaction between liquid 
subphase and solid substrate takes the form [3lj |: 



A 



12nd 2 



(37) 



Here A is the Hamaker constant and has the dimension 
of energy. It can be positive or negative depending on the 
properties of liquid subphase and solid substrate. If A is 
positive, van der Waals interaction leads to an effective 
repulsion between liquid-air and liquid-substrate inter- 
faces and favors a thicker liquid film, i.e. larger d. On 
the other hand, if A < 0, the liquid film can be unstable. 
Spontaneous fluctuations may rupture the liquid film via 
spinodal dewetting [40| . The value of A is typically in 
the range of 10" 20 J to 1CT 19 J From equation ([IS]), 
we have the critical wavenumber q c and wavelength A c : 
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(38) 



Thus q c and A c exist only in the case A > 0, i.e. for stable 
liquid subphase. The second-order buckling transition 
requirement (|36[) reduces to: 



A < 200 irB , 



(39) 



which doesn't depend on the thickness of the liquid sub- 
phase do- This independence of thickness can be under- 
stood by noticing that as we increase the thickness of 



Xr 



27rdo/V10 « 2.0 d . 



(40) 



The wavelength A c increases only gradually from the 
minimum because of the weak dependence on B/A in 
equation ([55)) . For practical purposes, the buckling 
wavelength is confined to scales of order do. In the 
microscopic range that we are discussing, do < lOOnm. 
Typically, the dimension of the Langmuir system L in 
horizontal direction is about 1mm, so the condition of a 
large enough system \q c L\ ^> 1 is satisfied. 



B. Charged Double-layer Interaction 

The charged double-layer interaction has the form [32j : 



4>(d) = ipo exp (— nd) 



(41) 



where 1/n is Debye screening length. The coefficient 4>o 
is a constant depending on zeta-potentials of two sur- 
faces and electrolyte concentration in between. For a 1:1 
electrolyte, 4>o can be written as: 



MkT lll2 C s 



(42) 



In equation (|42|) . T is the temperature and C s is the con- 
centration of electrolyte in the bulk. Moreover 7, (1=1,2) 
is related to zeta-potentials £oi of the two surfaces [32| : 



7l = tanh(eCo*/4fcT) . 



(43) 



The potential <fio is positive if two electrical surfaces 
have charges of the same sign and repel each other. In 
the case 4>o is negative, two electrical surfaces attract each 
other; this may rupture the liquid film via unstable modes 
of undulation. The above equation is valid in the weak 
overlap approximation , in which the overlap between 
electrical double layers is small (33|: 



k ■ d > 1. 



(44) 
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From equation (18]), we have the critical wavenumber q c 
and wavelength A c : 
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exp (—ndo) 
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(45) 



Thus A c and q c exist only in the case 00 > 0, i.e. for 
stable liquid subphase. The second-order buckling tran- 
sition condition (155)) takes the form: 



dn > — In 



K \K 2 B 

The smallest A c compatible with condition (|4"6"]l is 



(46) 



A c 



(47) 



It is on the order of the Debye screening length k . 
The wavelength A c is an increasing function of do. For 
example, we consider a negatively-charged mica sub- 
strate and a negatively-charged conventional liposome 
(lecithin/cholestrol 6:4 molar ratio). As an electrolyte, 
we consider NaCl at concentration C s — 0.001 mol/L. 
In this case, the Debye screening length is kT 1 ss 10 nm 
[32| . At pH 5.8, the zeta potentials of mica surface (34| 
and conventional liposomes are [351 ] : 



Comic = -104mV and Co I* = -20mV. (48) 
The value of cf>o is calculated as: 



0o « 2.2 x 10~ 4 (J/m 2 ). 



(49) 



By the estimation of the bending modulus B > 10/cT 
[361 ] . the condition (|4G[) reduces to: 



do > 8.1nm . 



Sharma Potential 



(50) 



The Sharma potential is used widely in studies of wet- 
ting phenomena between different liquid thin films and 
solid substrates [13, [M, [39[, e.g. water on mica. It in- 
cludes both the apolar(Lifshitz-van der Waals) and polar 
interactions and has the form: 



polar liquid and d c is the Born repulsion cutoff length 
[371 ]. In this case, the critical wavenumber q c and wave- 
length A c have the form: 



1 



6S AP d 2 , 6 d 



d B^ 
A (; = 2^d B 1 / 4 ( &S Ap dl 



P 

2—0. exp[{dc - do )/l]\ (52) 



exp [(d c - do) /I] 



1/4 



The existence requirement of second-order buckling 
transition yields: 



120 S AP d A c F > -S^exp 



doS p exp 



do 



d c - do 



2l/VB[6S Ap d 2 c l 2 



I 



.3/2 



(53) 



For example, we consider a lipid monolayer, with bend- 
ing modulus B > lOfcT, sitting on top of water with a 
mica substrate below. The numerical values of coeffi- 
cients are S AP = 20mN/m, S p = 48mN/m, I = 0.6nm 
and d c = 0.158nm [3a ]. The wavelength A c is an in- 
creasing function of d - In the microscopic range, where 
the thickness of water subphase do < lOOnrn, we get an 
upper bound for A c : 



A c ^$ 1-2 \im . 



(54) 



Obviously, the condition \q c L\ ^> 1 is guaranteed. Fur- 
thermore, the inequality ([53]) is true for any positive value 
of do- In other words, the buckling transition will al- 
ways be second-order if the sharma potential correctly 
describes the interaction between water subphase and 
a mica substrate. The following dimensionless quantity 
changes very slowly with the value of do: 



V 



B 1 / 4 QS Ap dl 



S^_df 

P 



exp [(d c - do) /I] 



-1/4 



In the range of do between lnm and lOOnm, r\ S [0.8, 
1.9]. Thus the critical wavelength A c can be written as: 



A c = 2irdoT] 
> 5.0 d . 



(56) 



For the estimation in the last step, rj takes the smallest 
value 0.8 in this range. 



m = S Ap .(^)+S p -exp 



d r — d 



I 



(51) 



where S AP and S p are apolar and polar contribution to 
the spreading coefficient. I is the correlation length for 



IV. NUMERICAL RESULTS 

In order to verify our results and explore the predicted 
wrinkling phenomena concretely, we have done a discrete 
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(a), flat chain of springs. (b). compressed chain of springs. 
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FIG. 2: Discrete model of one dimensional chain of springs 
used in the simulation, (a). Flat state with no compression, 
(b). Buckled state with large enough compression. 



numerical implementation of the system. The numer- 
ical simulation was carried out using the Mathematica 
program. We modeled the surfactant layer by a one- 
dimensional chain of nodes connected by springs with 
un-stretched length a and spring constant k. The un- 
stretched length a is set to be 1 in the simulation for 
convenience. In order to impose the inextensibility con- 
straint, we set the spring constant k to be a very large 
value. A bending energy of B9f ^ i+1 ^/2 is assigned to 
every pair of adjacent springs, where 6^(.j +1 ) is the an- 
gle between these two springs along the chain direction. 
The total bending energy is a sum of all pairs along the 
chain. In the numerical simulation we set the free liquid- 
air surface tension in our theory 70 = 0. The substrate 
potential is discretized correspondingly by replacing the 
integral with a summation along the chain. With no 
compression the chain adopts a flat configuration lying 
on the x axis. In the simulation, the first node is fixed at 
the origin. In order to reduce the influence of boundary 
effects in a finite system used in the simulation, we fix 
the z-coordinate of the last node to be zero, while its x- 
coordinate was determined by the amount of compression 
A, as shown in Fig[21 All the other nodes are movable 
both in x and z directions in the process of minimization. 
The total free energy g* of this discrete model takes the 
form: 

9* = n ex L + fb + ft + fk 

= \\ ex L Q +gl, (57) 

where fb, fi and fk are bending energy, substrate po- 
tential energy and elastic energy stored in the springs. 
The quantity g\ is the sum of these. For the flat refer- 
ence state, the above three energy terms are zero. So the 



FIG. 3: Configurations with 41 nodes as compression A 
changes. The box in the graph shows the amount of com- 
pression. Van der Waals interaction is used in this simula- 
tion. Parameter values: A = 200, B = 1, do = 3.78036, 
k — 10 11 . The predicted wavelength of wrinkling: A c = 10. 
The nonuniformity of the amplitude with 1% compression is 
due to nonlinear effects as approaching a possible wrinkle-to- 
fold transition The position of a fold could change for 
different initial conditions of minimization. 



total free energy for the flat state is: 

g*=Tl ex L 1 (58) 

where L is the total length of the chain. The compression 
is A = L — Lq. At each fixed amount of compression 
A, we minimized the value of <?i(A) and computed the 
smallest external pressure U m i n (A) needed to reach this 
compression via the following equation: 

n mm (A) = g ^ (A) , (59) 

We chose a value of do such that the predicted A c was 
commensurate with the system, namely A c = 10. Then 
starting from a finite amount of compression, e.g. A = 
1%, we gradually lowered the value of A. If our theory is 
correct, the chain will approach a sinusoidal shape with 
wavelength A c . Moreover II mi „(A) should approach II C 
as A goes to zero. The order of buckling transition can 
be deduced from the functional shape of II m i„(A). For 
a second-order or continuous transition, H m i n (A) is a 
monotonic increasing function of A. There is no jump 
when crossing the transition point. In the case of a first- 
order transition, n m j n (A) is not monotonic. It has a 
minimal value n^ in at a nonzero compression A*. As 
a result, as soon as H ex exceeds n^ in , the configuration 
will jump to this finite amount of compression showing 
the property of a first-order buckling. 
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FIG. 4: Configurations with 41 nodes as compression A 
changes. The box in the graph shows the amount of com- 
pression. Parameters values are the same as in Fig. 3. 



FIG. 5: Theoretically predicted IT C and II m i„ from simulation 
as functions of ratio A/B. B is fixed at 1. Solid line represents 
prediction from theory. Discrete points are values of II m i n 
from the simulation. 



A typical sequence of chain configurations with 41 
nodes as we changed the amount of compression is shown 
in Figj3] and [4] Figj5] shows the agreement between pre- 
dicted values of II C and Tl m i n from simulation for several 
cases. It is noticed that II C is a constant function of A, 
if the buckling wavelength A c is fixed. 

Verifying the predicted transition from continuous to 
discontinuous wrinkling proved to be more subtle. Even 
though the ratio A/B is nearly a factor of 2 above the 
predicted threshold, we didn't see any evidence of dis- 
continuous buckling using discrete model. To understand 
this required a second numerical method. As we discuss 
below, it reveals that the discontinuity is too weak to 
have been seen in the discrete model. 

The second numerical method is based on a continuous 
model. To simplify the calculation, we took a single 
sinusoidal mode as Ansatz: /isin(gx). For convenience, 
only one period of wrinkling is included in the continuous 
model, while a real system is composed of many copies of 
it. The wavelength A changes as we change the amount 
of compression A: A — X c — A. Using a similar approach 
as the discrete model described above, we can compute 
n m j„(A) for each amount of compression. The order of 
transition is still determined by the functional shape of 
H m i n (A). This time not only the values of II C , but also 
the order of buckling are in good agreement with our 
theoretical prediction. The paradox with discrete model 
is also explained. As it is shown in Figj6j the first-order 
transition is very weak in the case of van der Waals 
interaction. With A/B — 1000, H m in has a minimal 
value at 1% compression with a 0.6% change in Tl m in- 
Thus the buckling transition is first-order. However, 
such a small change in H m in cannot be detected in the 
above discrete model. Based on the above numerical 
results, our theory makes good prediction for the 



substrate-induced buckling transition. 



V. DISCUSSION 

In the previous sections, a substrate-bending model 
was constructed. Here we discuss the implications of 
this model. Using our model, we made prediction about 
wrinkling wavelength A c with large enough compression 
Tl ex > n c and the order of buckling transition. In order 
for our mechanism to account for wavelengths of hun- 
dreds of nanometers, the trapped fluid layer itself would 
have to be many nanometers thick. During the transfer 
of a monolayer to substrate, such thick fluid layers usu- 
ally exist. The compressive stress required to buckle the 
surfactant layer could be developed during this transfer 
and rapid drying after deposition. 

The wrinkling mechanism predicted here is expected 
in any supported monolayer or bilayer systems with suf- 
ficient compression. A well-defined wrinkling wavelength 
A c is given in terms of the thickness of the subphase 
do, the bending modulus of the surfactant layer B and 
functional form of the substrate potential 4>(d). Infor- 
mation about these microscopic variables is embedded 
in the experimentally observed wavelength. It is espe- 
cially useful if one can control the external pressure Tl ex . 
The properties of gravitational wrinkles have been ex- 
perimentally studied [9j . Gravitational buckling appears 
to give rise to a strongly first-order wrinkling-to-folding 
transition, which contrasts with our very weak first-order 
transition in the case of van der Waals interaction. Some 
other forms of substrate potential 4>{d) may give rise to 
a stronger first-order transition. 
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factant layer and etc. It will be of interest to see how 
such buckling occurs in time, and what counterparts of 
the wrinkling-to- folding transition might exist. 
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FIG. 6: Umin versus Log(A) for very weak first-order tran- 
sition in the case of van der Waals interaction. A — 1000, 
B = 1, do = 5.65295. The critical wavelength is A c = 10. 



These surfactant layers are potentially subject to an- 
other kind of instability different from the extensive wrin- 
kling investigated here. The boundary conditions may be 
such that the boundary region buckles while the bulk of 
the layer is still in a stable state. These boundary in- 
duced deformations of surfactant layers have been stud- 
ied [52, Ell . Such boundary buckling was an important 
factor in our discrete simulation. It prevented us from 
studying arbitrary wavelengths. Also our methods only 
allowed us to study the region of incipient instability. 
There may be interesting phenomena analogous to the 
gravitational wrinkling-to- folding transition that we have 
missed. 

Inextensibility of the surfactant layer is an important 
assumption in our theory. If this constraint were released, 
the system would have a compression mode as an extra 
degree of freedom to store elastic energy besides bend- 
ing mode studied above. As shown in the appendix, in 
the limit of small deformation approximation, finite com- 
pressibility influences only fourth order and higher terms 
in the free energy, so it doesn't change the expressions of 
the threshold external pressure II C and the critical wave- 
length A c . Moreover, if the system were not too far away 
from the transition threshold between first-order buck- 
ling and second-order buckling, the inextensibility would 
be a good approximation in the experiments of interest 
here. 



VI. CONCLUSION 

As supported monolayers and bilayers become more com- 
monly studied, we expect that the type of wrinkling pre- 
dicted here will be observed and used to infer local prop- 
erties, such as substrate depth, bending modulus of sur- 



APPENDIX A 

To estimate the influence of compressibility, we assume 
that t is the thickness of the surfactant layer. The bend- 
ing modulus B varies as t 3 , while the compressibility 
modulus K is proportional to t [44] : B/K ~ t 2 . The 
fourth order term of bending free energy takes the form: 



h 5B 

04 = — 



(Al) 



In the limit of small deformation approximation, 
derivatives of £ can be approximated as: 



£ ~ h/X c 



(A2) 



Inserting equation (|A2[) into the expression of 04, we get: 



hi 



5Bh 4 L 



(A3) 



The compression free energy takes the form: 



E, 



K 



V 



1/2 



dx . 



(A4) 



where K is the compressibility modulus. 77 is the percent- 
age change of the length of surfactant layer. Choosing the 
configuration of the surfactant layer just before buckling 
as a reference state, we have: 



n ea ^ lie 

K 



(A5) 
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where L c is the length of the surfactant layer just be- 
fore buckling transition. If buckling is not allowed, the 
compressive strain r\ is evidently given by the fractional 
decrease in L, viz (L c — Lq)/L c . However, if buckling 
occurs, this strain can only decrease. As a result, the 
following relation holds: 



V < 



(A6) 



where h is the slope amplitude. Thus, the expansion of 
compression free energy in terms of h, being quadratic in 
77, has only a fourth order term or higher. As a result, it 
doesn't affect the expressions of II C and A c . To the lowest 
order approximation, Ek can be written as: 



E k ~ c 4 



K V 2 L 



(A7) 



Comparing 64 and C4, we get a criterion of incxtensi- 
bility: 



(A8) 



where we have used the relation B /K ~ t 2 . By equation 
(|A5|) . we have: 



K 



n, 



n„ a 2 



K 



A 2 

h 2 -A 2 /K, 



(A9) 



where the expression of A 2 was given in equation (|29p . 
Inserting into the above criterion of inextensibility, we 
have: 



A 2 t 
K«X C 



(A10) 



As n ex approaches II C from above, A 2 can be approx- 
imated as: 



l - (B^(d )/^ 2 \d ) - 2{B^ 2 \d,)f/ 2 ) 
j l [B^(d )/cf > ^(do)~87r 2 B/X 2 c ) , (All) 



where we have used the expression of A c in the last step. 
We require A 2 > in order for the transition to be 
second-order. Thus, the first term in equation (|A11[) 
must dominate the second. However, if the system were 
not too far away from the transition threshold between 
first-order buckling and second-order buckling, two terms 
in equation (|A1 1|) would have comparable order of magni- 
tude. In such a case, we can simplify the criterion (|A10[) 
as: 



K 



B 

~xTk 



t 



i.e. 



A, 



< 1. 



(A12) 



The above criterion is always satisfied in the experimental 
systems that we are interested in. For example, suppose 
the monolayer thickness t is about 2nm and the wrinkling 
wavelength A c is more than lOOnm. In such a case the 
approximation of inextensibility is valid. In cases of very 
anharmonic potentials where 4>( 4 \do) ^> (j)^(do)/X 2 , the 
two terms in A 2 would not be comparable. Then the 
effects of compressibility could become significant. 
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